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Within the mode-coupling theory for ideal glass transitions, an analysis for the correlation func- 
tions of glass-forming systems for states near higher-order glass-transition singularities is presented. 
It is shown that the solutions of the equations of motion can be asymptotically expanded in polyno- 
mials of the logarithm of time t. In leading order, an ln(t)-law is obtained, and the leading corrections 
are given by a fourth-order polynomial. The correlators interpolate between three scenarios. First, 
there are planes in parameter space where the dominant corrections to the ln(t)-law vanish, so that 
the logarithmic decay governs the structural relaxation process. Second, the dynamics due to the 
higher-order singularity can describe the initial and intermediate part of the a-process thereby re- 
ducing the range of validity of von Schweidler's law and leading to strong a-relaxation stretching. 
Third, the ln(t)-law can replace the critical decay law of the /3-process leading to a particularly 
large crossover interval between the end of the transient and the beginning of the a-process. This 
may lead to susceptibility spectra below the band of microscopic excitations exhibiting two peaks. 
Typical results of the theory are demonstrated for models dealing with one and two correlation 
functions. 



PACS numbers: 61.20.Lc, 82.70.Dd, 64.70.Pf 



I. INTRODUCTION 



Within the mode-coupling theory for ideal glass tran- 
sitions (MCT), the dynamics of an amorphous system of 
strongly interacting spherical particles is described by M 
functions of time t, <j> q {t), q — 1,2, ... , M . These are au- 
tocorrelation functions of density fluctuations with wave- 
vector modulus q chosen from a grid of M values. The 
theory is based on a closed set of coupled nonlinear equa- 
tions of motion for the <j> q {t). The coupling coefficients 
in these equations are given in terms of the equilibrium 
structure functions. The latter are assumed to be known 
smooth functions of the control parameters of the system 
like, e.g., the packing fraction ip [[[]. The MCT equations 
exhibit fold bifurcations |2j at certain critical values of 
the control parameters, say at ip = (p c , describing a tran- 
sition from ergodic liquid dynamics for ip < tp c to non- 
crgodic glass dynamics for tp ^ ip c . The transition is ac- 
companied by the evolution of a slow stretched dynamics 
that was suggested as the explanation of structural re- 
laxation observed in glass-forming liquids. The leading- 
order asymptotic solutions of the equations for param- 
eters approaching the transition provide predictions for 
the universal properties of glassy dynamics j|. These 
predictions have been tested extensively against exper- 
imental data and molecular-dynamics-simulation results 
0, E|. The outcome of these tests qualifies MCT as a 
candidate for a theory of structural relaxation in glass- 
forming systems. 

It was shown that schematic MCT models exhibit also 
higher-order-bifurcation singularities like the cusp and 
swallow-tail bifurcation. The accompanying dynamics is 
utterly different from the one for the fold bifurcation. For 
example, in certain parameter regions, the leading order 
result reads (<f>(t) — /) oc — ln(i/r) ||. This logarithmic 
decay is equivalent to a susceptibility spectrum that is 



independent of frequency u>, "x!'(ui) x w j or to a 1/ f- 
noise fluctuation spectrum. There are corrections to this 
leading-order result which alter qualitatively the straight 
</>(i)-versus-logt lines or the plateaus of the x"(w)-versus- 
log(cj) plots. One needs to understand these corrections if 
one intends to get an overview of the relaxation scenarios 
for parameters near the higher-order-bifurcation points. 
It is he goal of this paper to provide such understanding 
by construction of a general theory for the logarithmic- 
relaxation law and its leading corrections. 

For parameters at a cusp singularity of schematic 
M = 1 models, the leading-order long-time decay fol- 
lows the law (4>(t) — f) oc l/ln 2 <. This law has been 
embedded in a leading-order description of the dynamics 
near the singularities in terms of multi-parameter scaling 
laws 0. It was shown by L. Sjogren that dielectric-loss 
spectra for certain polymers could be interpreted by this 
scaling-law description Q , and further work extended his 
analysis || [l0|, [H]. However, it was also demonstrated 
that the cited decay laws at the critical points have to be 
complemented by their leading corrections in order to de- 
scribe the numerical solutions of MCT equations within 
a time regime relevant for data analysis [7], [l2|. But, so 
far it was not possible to evaluate the corrections for the 
mentioned scaling laws. The results of this paper will be 
obtained along a different route of asymptotic expansion 
of the MCT solutions than followed in Ref. @. 

Logarithmic decay of correlations for glassy sys- 
tems has been observed, for example, in Monte-Carlo- 
simulation results for a spin-glass model |D| , for photon- 
correlation data from a dense colloidal suspension [fTi) , 
and for optical Kerr-effect data for a van der Waals liq- 
uid But the present work is motivated by three 
recent discoveries. First, density correlators 4> q {t) mea- 
sured by photon-correlation spectroscopy for colloids of 
micellar particles demonstrated logarithmic decay within 
time windows of two orders of magnitude in size p6[. 



2 



Second, the MCT equations for a system whose struc- 
ture was described by Baxter's model for sticky hard 
spheres exhibit cusp bifurcations [Jl?], These find- 

ings have been corroborated by a comprehensive analy- 
sis of the glass transitions of a square- well system [[l9| . 
Third, logarithmic decay extending over three decades 
in time was found in a molecular-dynamics simulation 
for a system with an interaction given by a strong re- 
pulsion complemented by a short-ranged attraction p0[ . 
One concludes that higher-order bifurcation singularities 
are not restricted to schematic models and that there are 
reasons to suggest the search for such singularities in col- 
loids with short-ranged attraction. It is the aim of this 
paper to provide a detailed discussion of the qualitative 
features, which are characteristic for the relaxation in 
systems near higher-order bifurcations. A set of general 
formulas shall be derived, which could be used as ba- 
sis for a quantitative analysis of future experiments and 
simulation studies. 

The paper is organized as follows. In Sec. II A, the 
known general MCT equatio ns fo r structural relaxation 
are formulated. Then (Sec. II B), these equations are 
rewritten in a form which is suited as a basis for an 
asymptotic solution near bifurcation singularities. Sec- 
tion H| presents the theory for the logarithmic relaxation 
for MCT models dealing with a single correlator and in 
Sec. ^ quantitative results are discussed for a cusp sin- 
gularity in an M = 1 model. Sectio n M presents the 
theory for the general case and in Sec. [VJfurther results 
are discussed for relaxation near a swall ow-tail singular- 
ity for an M = 2 model. Section VII summarizes the 
findings. 



II. BASIC EQUATIONS 

A. The equations for structural relaxation 

MCT is based on two sets of equations. The first one 
consists of the exact equations of motion for the M den- 
sity correlators 4> q (t), q = 1,2, . . . , M, derived within the 
Zwanzig-Mori formalism 

d?<p q (t) + Sl%{t) + f M q (t - t')d t ,<p q (t')dt' = . 
Jo 

(la) 

The initial conditions read <p q (t = 0) = 1, dt<fi q (t = 0) = 
0. The positive tt q are characteristic frequencies and the 
kernels M q (t) are fluctuating- force correlators plf . In 
colloidal suspensions, there are contributions to the force 
due to interactions of a colloid particle with the solvent 
particles. These fluctuate on a time scale much shorter 
than the one relevant for the motion of the mesoscopic 
colloid particles. Therefore, one can approximate the cor- 
responding contributions to the kernel by a white-noise 
term, and this leads to a friction force v q dt<j) q {t), v q > 0. 
Compared to this term, one can neglect the inertia term 



dt<fi q (t). One arrives at the analogue of Eq. ( |Ta| ) for col- 
loids, i.e. at an equation of motion where the underlying 
dynamics is Brownian rather than Newtonian: 



v q d t <i>S) + + / M S - t')d t ^ q {t')dt' = o . 

Jo 

(lb) 

The initial conditions are cf> q (t = 0) = 1. The ker- 
nels are split into regular ones, M q es (t), and so-called 
mode-coupling kernels m q {t) describing the cage effect: 
M q (t) = M^ cg (<) + tt q m q (t). The regular terms describe 
normal liquid effects like binary collisions in conventional 
liquids or hydrodynamic interactions in colloids. The 
crucial step in the derivation of MCT is the applica- 
tion of Kawasaki's factorization approximation in order 
to express m q {t) as functionals T q of the M correlators 
4>k(i) called mode-coupling functionals. For simple sys- 
tems, they are quadratic polynomials, whose coefficients 
are given in terms of the equilibrium structure functions 
El p2). The functionals depend smoothly on, say, N con- 
trol parameters to be combined to a control parameter 
vector V = (V\,... ,Vn)- In conventional liquids, the 
packing fraction and the temperature may be the control 
parameters. In colloids, one of the control parameters 
may be the attraction strength, which can be changed 
by modifying the solvent. Let us note the functionals by 
FgWJk]. For < \fk\< 1, k = 1,2,... ,M, they can 
be written as Taylor series with non-negative coefficients. 
Thus, the second set of MCT equations is 



m q {t)=F q [V,(t>k{t)] 



(2) 



Specifying the functionals T q , the regular kernels 
M q cg (t), and the frequencies h q or v q , Eqs. ( pi] ) or 
Eqs. (lb) together with Eq. (||) are closed. In the present 



paper, a topologically stable singularity of Eqs. ([!]) and 
(^) shall be discussed in full generality. Therefore, mi- 
croscopic details are not of concern. 

It will be convenient to discuss dynamics in the do- 
main of complex frequencies z, Imz > 0. This can be 
achieved by Laplace transformation of functions of time, 
say F(t), to functions of z denoted by C\F(t)](z) = 
i J °° e~xjp(izt)F(t)dt. Equations ( p"aj ) and (|) together 
with the initial conditions are equivalent to a fraction rep- 
resentation of 4> q {z) = C[4> q (t)](z) in terms of M™ g (z) = 
C[M q ^(t)]{z) and m q {z) = C[m q (t)}(z), 

4> q {z) = -l/{z - Q 2 J [z + M^(z) + n 2 q m q (z)] } . (3a) 

The analogue for the colloid dynamics is derived from 
Eq. ph: 



<P q {z) = -l/{z - [iu q + M™z(z) + n 2 q m q (z)]} . 

(3b) 

There are two possibilities for the solutions of the pre- 
ceding equations. Either, all long-time limits of the cor- 
relators vanish as expected for an ergodic system. States 
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with such control parameters 7 are referred to as liq- 
uids. Or, there may be nonvanishing long-time limits 
fq — — * < f q < 1, as expected for non-ergodic 
systems. States with such control parameters 7 are re- 
ferred to as glasses, and f q is called the glass-form fac- 
tor. Changing 7, it may happen that there are values V c 
where one changes from a liquid to a glass — these are 
the ideal liquid-glass transitions discussed within MCT 
0. For glass states, 4> q (z) exhibits a zero- frequency pole, 
<j) q (z —>■ 0) ~ —f g /z. Because of Eq. (||), a similar state- 
ment holds for the kernel, m q (z — > 0) ~ —J- q [V,fk]/z. 
Hence, for frequencies tending to zero, kernel m q (z) be- 
comes arbitrarily large compared to the term z + M q cg (z) 
or the term iv q + M q es (z), respectively. Because of con- 
tinuity, also for states with control parameters near the 
ones for glass states Q q m q (z) is very large compared to 
z + M r q C6 (z) or to iv q + M* es (z). Under the specified 
conditions, Eqs. (pa, b) simplify to S 



ct> q {z) = -l/[z-l/m q {z)] 



(4) 



This equation exhibits most clearly the fraction represen- 
tation of correlators, which is the essence of the Zwanzig- 
Mori theory. It shows that there is a one-to-one corre- 
spondence between density fluctuations and force fluctu- 
ations. There is no separation between the time scales 
for the particle motion within cages and for the particles 
forming the cage. Therefore, the correlators <fi q (t) and 
the kernels m q (t) have to be calculated self-consistently 
- this is the essence of MCT. Equation (||) is scale in- 
variant. With 4> q (t), also (j) q (t) — <f) q {x-t) is a solution for 
any x > 0. The scale for the high frequency dynamics 
is determined by the transient motion and this is gov- 
erned by M q cg (t) and £l q or v q . But these quantities do 
not occur anywhere in Eq. (|J). Thus, Eq. (Q) can fix the 
solution only up to some time scale. 

Equations (|J) and (^) are the MCT equations for struc- 
tural relaxation. In particular, they are the basis of the 
asymptotic expansions for the long-time dynamics for 
control parameters near bifurcation points. 



B. The equations for structural relaxation near 
glass-transition singularities 

In this section, the concept of a glass-transition singu- 
larity shall be reviewed. The equations of motion shall 
be rewritten in a form where the small quantities, which 
characterize the relaxation near such a singularity, ap- 
pear transparently. 

To simplify the notation of the following equations, the 
Laplace transform shall be modified by a factor of {—z)\ 



S [F{t)\ (z) = {-iz) / exp(izt) F(t) dt . (5) 
Jo 

This linear mapping of functions of time to functions of 
frequency leaves constants invariant: F(t) — c implies 



S[c](z) = c. Let F(t) = (X*(t)X)/(k B T) denote a cor- 
relation function for variable X determined for tempera- 
ture T. Then, the dynamical susceptibility for frequency 
lu can be written as x{oj) — F(t = 0) — S [F(t)j {w + iO) 
pl[ . Thus, S [F(t)] (z) denotes the non-trivial part of a 
dynamical susceptibility. 

Equation ([IJ) can be rewritten as 

s iMt)] (*)/{i - s IM*)] (*)> = s fa [V, Mm (*) ■ 

(6) 

Equation <p q (t — + oo) = f q is equivalent to S [</> q (t)] (z — > 
0) = f q . Similarly, one obtains S [FqjV, (f>k(t)]] {z — > 0) = 
F q [V, fk}- The z -> limit of Eq. @ yields a set of M 
implicit equations for the M glass- form factors f q |jj : 

(7) 



h/^-f q ) = F q [v,f k ] 



This equation may have other solutions, say f q . The 
glass-form factor is distinguished by the maximum theo- 
rem: /, < /„ q = 1, . . . , M ||. 

Let V c denote some reference state. The long-time 
limits of the correlators for this state shall be denoted 
by f q . The correlators shall be written in terms of new 

functions <j> q {t) 



4>q{t)=f C q + {l-fl)^S) ■ 



(8) 



The functional J- q [V, <fik(t)] can be rewritten as a Taylor 
series in <pk(t), using the coefficients 



.4 



(n) 
qki- 



XV) 



i 

n 

x (1 



^i~f c q ){d n ^ q [v,m /dr kl ---dr k j 
-fkj---^-fkj- ( Qa ) 



These shall be split in the values for the reference state, 



A 



(n)c 
qk\ ---kn 



.(») 

qk± •••k. 



(V c ), and remainders: 



A 



U n (V) = A 



(n)c 
qk\ ■ --kn 



A [n) 

^■qki—kr. 



(V) . (9b) 



Let us consider a path in control-parameter space given 
by V{e) = (Vi(e), . . . V N {e)). The N components of 7(e) 
are smooth functions of the path parameter e, and the 
tangent vector dV(e)/de must not vanish. Let us choose 
7(e = 0) = V c , so that e can be considered as a distance 
parameter specifying the neighborhood of V c . One gets 
7(e) = Te + 0(e 2 ), with T = dV{0)/de being the tangent 
vector of the path at 7 C . The mode-coupling functional 
is a smooth function of 7, i.e. 



.4 



(n) 
qki . 



, (7)= 0(e) 



(9c) 



The details of the e-dependence of the coefficients are 
not important. The parameter e is introduced mainly as 
a means for bookkeeping in the following expansions in 
7— 7 C . Expanding the left hand side of Eq. (||) in powers 

of S (j) q (z), one can rewrite this equation in the form 



Sqk - A 



(l)c 
qk 



S[Mt)}(z) = Jq(z) 



(10a) 
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j q (z) = AW(v) + AU(v)S[4> k (t)](z) 

oo 
?i=2 

-S[$ q (t)] n (z)} ■ (10b) 

Here and in the following, summation over pairs of equal 
labels k is implied. These equations are equivalent to 
Eqs. (Q) and ([I|) for the structural relaxation. The 
small quantities to be used for the asymptotic solution 
are the coefficients k and the functions 4> q (t) or 

S (j) q (t) (z), respectively. 

Specializing Eqs. jl0| ) to the i 
equation for f q — 4> q {t — > oo): 



calculating M correlators to the one of calculating the 
projection of the correlators on the dangerous eigenvec- 
tor of the above defined Jacobian. Second, one has to 
solve the equation for the projection. In this section, the 
second problem will be studied, which is equivalent to a 
discussion of M = 1 models. 



A. Classification of glass-transition singularities 

One component models deal with a single correlator 
<j>(t), a single glass-form factor /, etc. All matrix indices 
can be dropped in the formulas of Sec. II. The one-by- 
limit, one gets the one matrix A^ c is identical with its maximum eigenvalue 



[S qk 



qk iJ*-A q Hv) + A<V(v)f k 

oo 



+ 



(11) 



This is a rewriting of Eq. (|t|) so that small deviations of 
f q from and V from V c are explicit. The M by M 

matrix \5 q k — A^? c ] is the Jacobian of the set of implicit 
equations (0) for the reference solution at V = V c . 
This Jacobian consists of a unit matrix 8 qk and a matrix 

^qk° °f positive elements. The Frobenius theorems im- 
ply that, generically, this matrix has a non degenerate 
maximum eigenvalue E c > 0. All other eigenvalues have 
a modulus smaller than E € j2^| . It is a subtle property of 
MCT that E c ^ 1 Q . If E c < 1, the implicit-function 
theorem guarantees that all states V for sufficiently small 
e are states whose long-time limits f q = f q +(l — f q )f q 
depend smoothly on e. In the following, the reference 
state V c shall be specialized so that 



E c 



1 



(12) 



In this case, V c is a bifurcation point of Eq. (j^). The f q 
are singular functions of e for e — * 0, and therefore V c 
is referred to as a glass-transition singularity. Since E c 
is non-degenerate, the possible bifurcations are from the 
so-called cuspoid family Ai, I = 2,3, ... . The bifurcation 
singularity A\ is topologically equivalent to that for the 
zeroes of a real polynomial of degree I Q . The A 2 , also 
called fold bifurcation, is the generic singularity obtained 
by varying a single control parameter. The liquid-glass 
transition of MCT is of this type. In this paper, the 
dynamics near a higher order singularity Ai, I ^ 3, shall 
be analyzed. 



III. RELAXATION DESCRIBED BY 
ONE-COMPONENT MODELS 

It will be shown in Sec. [v| that each iteration step of 
the asymptotic solution of the equations of motion splits 
into two parts. First, one has to reduce the problem of 



E c . Because of Eq. (JIJ), the left-hand side of Eq. (|TJ 
vanishes. The equation for / reads ei(V) + C2(V)f 
^2 n>2 [A^ {V) — /in]/" = 0, where the abbreviations are 



i>2 

used 



Mn = l-A (n)c , e n (V) = A^(V) , n=l,2, 



(13) 



The singularity exhibited by / for e tending to zero de- 
pends on the number of successive vanishing coefficients 
fi n . A singularity of index 1,1^2, shall be defined by 



Mi = M2 



fii-i = o , m ^ o 



(14) 



The equation for / reads 

W/' =ei-i(V)f l ~ 2 + e^ 2 (V)f 1 - 3 + ■■ ■ + ei (V) 
+{ei(V)f l - 1 +e l+1 (V)f l + [A (n) {V) 



- 1 



n 

(15) 



The implicit-function theorem can be used to show that 
there is a smooth invertible transformation of the I 
variables ^ei, e 2 , ■ ■ ■ , which eliminates the curly 

bracket in Eq. (|l~5|). Thus, the singularities described by 
this equation are topologically equivalent to the ones de- 
scribed by the first line, i.e. by the zeros of a polynomial 
of degree I. In Arnol'd's terminolog y p[ , such singularity 
is referred to as A;. Because of Eq. (pcj), the e n (V) are of 
order e and shall be referred to as separation parameters. 

The simplest glass-transition singularity is the A 2 . In 
this case, there is only one relevant control parameter 
ei(V). One infers from Eq. ( |l5|) that there is a discontin- 
uous change of / at the surface specified by e\{V) = 0. 
The bifurcation dynamics is characterized by power-law 
decay and there appear power-law dependencies of the 
relaxation scales on |ei(V)|. All exponents in these laws 
are to be calculated from A = 1 — fi 2 , which is called the 
exponent parameter p2] . The transition surface has a 
boundary that is determined by A = 1, i.e. by [i 2 = 0. 
These endpoints are the higher-order singularities. The 
A3 and A4 are also referred to as cusp and swallow-tail 
singularities, respectively. 
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B. Equations for an asymptotic solution 



Let us specialize Eqs. (10a) and (10b) for M = 1. Let 
us also express the coefficients A^ n> (V) in terms of fj, n 
and e n (V): 



= ei00 
+e 2 (V)Sm)]{z) 



+ 



(1 
(1 
(1 



M3 )5[0 3 (t)](z) 

» A )s[i 4 (t)](z) 



(16) 



This suggests an expansion of the solution in powers of 
lei 1 / 2 . With G^(t) = 0(|e|™/ 2 ), let us write 



4>(t) = G (1) (i) + G (2) (i) + G (3) (i) 



(17) 



The first line of Eq. ( jig ) is of order |e| and it provides 
a nonlinear integral equation for G^(t). The contribu- 
tions to this line which are of order |e| 3//2 together with 
the leading terms of the second line provide a linear inte- 
gral equation for G^ 2 '(t), etc. This procedure will yield 
the desired asymptotic expansion provided the indicated 
integral equations define meaningful solutions. This is in- 
deed the case as shall be demonstrated below by explicit 
construction of the G^ n \t). To proceed, the following 
discussion shall be restricted to higher-order singularities 
by requiring 



M2 = 



C. The leading contribution 



The equation for the leading contribution to the cor- 
relator at a glass transition Ai with I ^ 3 reads 



(V) + S[G^\t)](z) - S[G W (i)] 2 (z) = . (19) 



The formulas for the Laplace transforms of the loga- 
rithm and its square imply <S[ln(t)](z) = \n{i/z) — 7 
and S[\n 2 (t)]{z) = ln 2 (i/z) - 2 1 \n{i/z) + 7 2 + (7r 2 /6), 
wher e 7 = 0.577 ... is Euler's constant (cf. Appendix [A|, 
Eq. Q). Hence, Eq. (|l|) is solved by -Bhx(t) if 
ei(V) + (B 2 ir 2 /6) = 0. Since the correlators are mono- 
tonically decreasing functions of t p4fl , one must require 
B > 0. One concludes that a solution is given by 



G« (t) = —B ln(t) , B = Vl-eei^/Tr 2 
provided the control parameters V obey 
ci(V) <0 . 



(20) 



(21) 



A more general solution is C?W(t) = C/W(t) + c, where 
c can be any real constant. Introducing x = exp(— c/B), 
one gets & x \t) = G^(xt). Thus, the generalization 
is the one implied by the scale invariance of the basic 
Eq. (Q). It shall not be considered here. Rather it shall 



be accounted for at the end of all calculations by rescaling 
t to t/r. Ignoring corrections of order |e|, one derives 
from Eqs. (||), ((ij^ and (^0|) the leading approximation 
for the correlator 



J>{t) =/ c -(l-/ c )Bln(t/r) 



(22) 



Let us anticipate that the smooth function ei(V) is 
generic for V near V c and has a nonvanishing gradient. 
Then, ei(V) = defines a smooth surface through V c in 
the control-parameter space. It separates the neighbor- 
hood of the glass-transition singularity V c into a strong- 
coupling side where £i(V) > and a weak-coupling side 
where £i(V) < 0. The results of this paper refer to the 
latter regime. 



D. The leading correction 

In order to solve Eq. (0) up to order |e| 3 / 2 , one 
has to incorporate from the first line the contribution 
2S[G« {t)G^ (t)] (z) - 2S[G«] (z) S[G (2 ^ (t)] (*), one has 
to evaluate the second line with <j> replaced by G^(t), 
and one can ignore all other terms. Hence, the equation 
for the leading correction G^ 2 '{i) can be written in the 
form 



(23) 



(18) Here, the linear integral operator T is defined by 



T [<?(*)] (z) = S {Ht)G(t)} (z) - S [ki(i)] (*)S [G(t)] (z) , 

(24) 



and the inhomogencity of Eq. ( |23| ) reads 

fW(z) = {e 2 (V)S[G^(t)](z) - fi 3 S[G^ 3 (t)](z) 

+2C{S[G^ 3 (t)}(z)~S[G^(t)} 3 (z)}}/(2B) . 

(25) 

A factor 2£ has been introduced for later convenience. 
For the study of M = 1 models, one has to substitute 
C - 1/2. 

The solution of Eq. fl9| ) was built on the equations 
T[c](z) = 0, T[\nt] \z) = tt 2 /6. These formulas are 
generalized in Appendix ^ by constructing polynomials 
p n (x) of degree n 1 with the properties: 



Pn(x) = b nt iX + b n , 2 X H \-b n ,n-lX n + X n , 



T[ Pn (Ht))} (z) = n(ir 2 /6) ln"" 1 ^) 



(26a) 



(26b) 



These polynomials are a convenient tool to solve the 
equation 



T[g(t)} (z) = f(z) 



(27a) 
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for inhomogeneities f(z) which are polynomials in 
hx(i/z), 



f(z) = ^2a n \n n (i/z) 



n=0 



Obviously, 



m+l 



g(t) = J2 K-l/(n7T 2 /6)] Pn(Ht)) 



(27b) 



(27c) 



Using Eq. (gOj) and applying Eqs. (|Al|) and ( |A2|) for 
the evaluation of the transformations of the p ower s of 
ln(i), one can write f^ 2 '(z) in the form of Eq. ( 27b ) for 
to = 3. The coefficients are linear functions of ei(V) and 

a = [(6CA 2 )(r 3 - if) - (3 W /7r 2 )r 3 ]ei(V) - (T 1 /2)e 2 (V) 

(28a) 

ax = [3C - (9M3A 2 )r 2 ]ei(y) - (l/2)e 2 (V) , (28b) 

a 2 = -(9/i 3 A 2 )r iei (T/) , a 3 = -(3/i 3 A 2 )ei(U) . 

(28c) 

Here, Tk = d k T(l)/dx k denotes the fcth derivative of the 
gamma function at u nity. One concludes that G^ 2 '(t) = 
g(t), where Eq. (27c) is to be used with to = 3: 



3=1 



(29a) 



The coefficients are derived with Eqs. (A7a-c): 

B x = (0.44425C- 0.065381^3)ei(^) - 0.22213e 2 (U) , 

(29b) 

B 2 = (0.91189C + 0.068713^ 3 )ei(V0 - 0.15198e 2 (U) , 

(29c) 

B 3 = -0.13504^ 3 ei (V) , Bi = -0.046197/^ 3 ei (V) . 

(29d) 

Dropping corrections of order |e| 3 ^ 2 , the solution up to 
next-to-leading order reads 

0(t) - f = (1 - f) [(— -B + B X ) ln(t/r) + B 2 ln 2 (t/r) 

+B 3 \n 3 (t/ T ) + B± ln 4 (t/r)](30) 

A singularity Ai with I ^ 4 implies /i 3 = 0. In this case, 
the formula simplifies because B 3 — B4 = 0. 

The described procedure can be continued. To solve 
Eq. (pit) up to order e 2 , one derives the analogue to 
Eq. @: T[G( 3 HiWz) = f {3) (z). Function f^(z) has 
the form of Eq. (E7q) with to = 6, where the coefficients 
<Zj depend on the parameters £i(U), e 2 (U), e 3 (V), /i 3 and 
/14. As a result, one gets 

G* 3 >(t)=£C;W (»/*), (31) 

3=1 

where Q = C(|e| 3/2 ). 



IV. RESULTS FOR A ONE-COMPONENT 
MODEL 

The simplest example for a generic cusp bifurcation is 
provided by an M = 1 model with the mode-coupling 
functional J-[V, f] — v%f + v 3 f 3 . This model was derived 
originally within a microscopic theory of spin-glass tran- 
sitions |2q] . It shall be used here in order to demonstrate 
several implications of our theory Let us use the model 
with a Brownian microscopic dynamics so that Eqs. ( |lb| ) 
and (H) specialize to 

T 1 d t 4>(t)+(t)(t)+ [ m{t-t')d t «j>(t')dt' = , (32a) 
Jo 

m(t) = VKt>(t) + v 3 (/) 3 (t) . (32b) 

' The two coupling constants v\ and V3 ^ are consid- 
ered as the components of the control-parameter vector 
V = {vi,v 3 ). 

Figure 1 reproduces the phase diagram fq, [22|. It 
is obtained from the largest of the solutions for f c of 
Eq. (§, i.e., v c J c + v c 3 f c3 = /7(1 - n, and Eq. ©, 
i.e., vl + 3wf/ c2 = 1/(1 - f c ) 2 , < f c < 1. There are 
two transition lines. The first one is the straight hori- 
zontal line of degenerate A 2 bifurcations: — 1, ^ 
v 3 ^ 4, f c — 0. Crossing this line by increasing vi, / = 
<p{t — > 00) increases continuously. The second one is the 
smooth curve of A 2 singularities V^ 2 '° shown as a heavy 
full line. It starts at v[ 2)c = 0, v (2)c = 27/4, f^ c = 2/3. 

. , f( 2 >° decreases along the line. For 
4 one gets f^ c = 1/2. Decreasing fW c 



(2) 

With decreasing ir ' 



,(2)c 



1, v. 



(2)c 



further, the line reaches the endpoint that is marked by 
a circle. This is the A 3 singularity V c specified by 



9/8, 



27/8, / c = l/3, /i 3 = l/3. (33) 



The two separation parameters are obtained from 
Eqs. (|9b|) and (|l^) as linear functions of the parameters 
differences Vx 3 = ^1,3 — 



ei(U) = (2/81)[90i + 1> 3 ], e 2 (V) = (4/27)[3fi x + fig] . 

(34) 

These formulas determine the coefficient B in Eq. ( |20| ) 
and B1-B4 in Eqs. (p9|). The scales r for the results 
in Eqs. (p2|), (|30|), and ([n]) are determined as the time 
where the correlator crosses the critical form factor: 

m = r. _ 

The dominant deviation of the correlators from the 
logarithmic-decay law, Eq. (|22]), is caused by the term 



B 2 In (t/r) in Eq. (p0[). Thus, the logarithmic-decay law 
is exhibited best for states V with B 2 = 0. This line is 
shown dashed-dotted in Fig. Figure || demonstrates 
the evolution of the dynamics upon shifting states on 
this line towards the A 3 singularity. The ln(t/r) interval, 
where Eq. (j|^) or (^) describe the correlators within an 
error margin of 5% is marked by closed or open symbols, 
respectively. For n ^ 2, these intervals increase with 
decreasing V — V c . 
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FIG. 1: Phase diagram for the one-component model defined 
in Sec. The horizontal light full line marks the liquid- 
glass-transition curve connected with a continuous variation 
of the glass-form factor. The heavy full line presents the set 
of A2 singularities, that ends at the A3 singularity marked 
by a circle. The dashed straight line describes the points of 
vanishing separation parameter ei and the dashed-dotted one 
the points of vanishing coefficients B2- Crosses with labels n 
and triangles with label n' denote states discussed in Figs, 
and Figs.]^, [| respectively. 




FIG. 2: Correlators <f)(t) for the one-component model defined 
in Sec. The states are located on the line B2 = with 
coupling constants: v{-vi = 0.9298/4", wg-v 3 = 3.3750/4", 
n = 1, . . . , 4, marked by crosses in Fig. ^. The full lines are the 
solutions of Eqs. (^2^,, b) with n = 1 as unit of time as also 
in the following figures. The dotted straight lines exhibit the 
leading approximation, Eq. the dashed lines the leading 
correction, Eq. (^o|). The filled and open symbols, respec- 
tively, mark the times where these approximations deviate 
from the solution by 5%. The dotted line marked by D is the 
Debye law exp[— t/n]. The horizontal line shows the critical 
form factor f c = 1/3. 



There are two peculiarities concerning the range of ap- 
plicability of the asymptotic expansions. First, it can 
happen that for sufficiently large e the range shrinks 
if one proceeds from the leading approximation to the 
next-to-leading one as is demonstrated in Fig. || for the 
n = 1,2 results. This is caused by a cancellation of er- 



rors due to neglecting the i?i-correction in the prefactor 
of the \n(t/r) term in Eq. (|30| ) and due to neglecting 
the terms proportional to B3 and B4. This peculiarity 
would disappear if the tolerated error margin were de- 
creased sufficiently below the 5% used. Second, for small 
V — V c , the interval of decay for (p(t) below the criti- 
cal form factor / c that is described by the asymptotic 
expansion shrinks with decreasing separation. This is in- 
ferred by comparing the n — 3 with the n — 4 results. 
The reason is the following. The correlator <p(t) decreases 
monotonically towards its long-time limit / |2J] . But the 
interval f c — f = —f shrinks for e — > 0, since Eq. ( |T5| ) 

implies -/ = (- ei (V0//i 3 ) 1/3 [l + 0(e 1/3 )]- 

Figure ^| demonstrates that the transient regime ex- 
tends to about t/n = 1. For vanishing mode-coupling 
functional, the correlator describes a Debye process: 
4>(t) = exp(— t/n)- Mode-coupling effects cause a slower 
decay for t/n ^ 1. But for V close to V c , the transient 
dynamics is rather insensitive to changes of the coupling 
constants. There is a crossover interval, say n < t < t*, 
before the decay of <fr(t) towards f c can be described by 
the ln(t/r) law. The beginning r* of the range of valid- 
ity of Eq. ( |2^ ) is indicated by the filled symbols. There 
are two subtleties demonstrated for n > 2. First, the 
time t* increases upon approaching the A3 singularity, 
and therefore the decay interval 4>(t*) — f c which is de- 
scribed by the logarithmic law shrinks with decreasing 
separation parameters. The control-parameter sensitive 
structural relaxation is governed by the two time scales 
t* and t. Both times become large, but r/r* becomes 
large as well for e — > 0. Second, the beginning of the 
range of applicability of the leading correction Eq. ( |30| ) 
is control-parameter insensitive, as is shown by the open 
symbols on the short-time part of the decay curves. 

The explanation of the findings in the preceding para- 
graph is based on the fact that for every fixed finite 
time interval, the MCT solutions are smooth functions 
of the control parameters [E4J. Therefore, for e tending 
to zero, the correlator for state V has to approach the 
correlator for state V c , the so called critical correlator. 
The latter is shown as the heavy full line with label c in 
Fig. |j| Thus, for every time interval ^ t i ma x and 
every error margin, there exists an e*, so that 4>(t) agrees 
with the critical correlator within the error margin for 
all |e| < e* and all ^ t ^ t max . This feature is demon- 
strated in Fig. H by the two curves with label n = 4 and 
e x ^ 0. They refer to states with v x -v{ = ±0.9298/4 4 
and U3 — U3 = ±3.3750/4 4 . These correlators are very 
close to the critical one for t ^ t max ; i max ~ 325 is in- 
dicated by open circles in Fig. ||. The critical correlator 
does not exhibit an ln(t/r)-part. Thus, the time r* for 
the onset of the description by Eq. ( pp| ) has to increase 
beyond any bound if e tends to zero. The asymptotic ex- 
pansion in Sec. Ill was based on \4>(t) — f c \ to be small. 
This condition is fulfilled for the critical correlator if the 
time is sufficiently large, since 4>(V c ,t) decreases mono- 
tonically to / c . Hence, for e — > there must appear the 
increasing time interval n < t < t* where the asymp- 
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FIG. 3: The heavy full curve with label c is the solution of 
Eqs. (j^a, b) at the As singularity, V = V c , called the crit- 
ical correlator. The thin full curves with labels n = 1 and 
4, ei < reproduce two of the solutions discussed in Fig. ^| 
The dashed lines reproduce from Fig. ^ the corresponding 
approximations by the next-to-leading order, Eq. (pOj). The 
dashed-dotted curves extend the asymptotic expansion by in- 
cluding the second order corrections, Eq. (|3l|). The hori- 
zontal line marks the critical form factor / =1/3. The 
full curve marked with n = 4, ei > refers to a state, 
vi-vf = 0.9298/4 4 , vs-v$ = 3.3750/4 4 . The circles mark 
the times where the correlators n — 4, e\ ^ deviate by 2% 
from the critical one. 



totic expansion describes the critical correlator. The e- 
dependence due to the separation parameters 6 n (V) and 
the e-dependence due to the time scale r cancel to pro- 
duce the critical correlator outside the transient regime 
and prior to the onset of the ln(t/r) law. This is shown 
most clearly by the dashed-dotted lines in Fig. |[ They 
exhibit the result of the asymptotic expansion up to the 
second correction given by Eq. (|3l]). They describe the 
complete decay for t/ri > 1 except for the final exponen- 
tial approach towards / = <fi(t — > oo); and this for states 
as far from the critical point as given by the one with 
label n=l (cf. Fig. [lj). 

Figure |] exhibits the dynamical susceptibilities x( w ) = 
1 - S[4>{t)](u; + iO) = x'M - «x"M for the states dis- 
cussed above. Without interaction effects, the suscepti- 
bility spectrum were a Debye peak x'b = Cdwtd/[1 + 
(wTfl) 2 ] with Cd — 1, To = t±. Such Lorentzian spec- 
trum is added to the upper panel as dotted line with 
label D, where Co and td are fitted to the maximum 
of the critical susceptibility spectrum. This shows, that 
the spectral peaks near w = 1, in particular their high- 
frequency wings, are due to the transient dynamics. How- 
ever, the low-frequency wings of the peaks are enhanced 
relative to the Debye spectrum and they are stretched to 
lower frequencies due to the critical relaxation within the 
interval 1/r* < lu < 1/t±. It is the structural relaxation 
near the ^3 singularity which causes the skewed shape of 
the spectral peaks. The leading approximation, Eq. (22), 
implies constant-loss plateaus: ~x"(ijj)/(l — f c ) = irB/2. 
However, this formula describes the plateau only for 
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FIG. 4: Susceptibility spectra x"(u;) and reactive parts of the 
dynamical susceptibility x'( w ) f° r the one-component model 
defined in Sec. IVL The full lines with labels n and c corre- 
spond to the correlators with the same labels shown in Figs. || 
and ^. The dotted straight lines, the dashed lines, and the 
dashed-dotted lines are the leading approximation, Eq. 
the leading correction, Eq. (3^), and the second correction, 
Eq. (|3ll), respectively. The filled, open, and half-filled symbols 
mark the frequencies where the corresponding approximation 
deviates from the spectrum by 5%. The dotted line with la- 
bel D exhibits a Debye spectrum Cd(ojtd) /[l + (lotd) 2 ] with 
Cd ~ 0.2503 and td/tl = 0.670 fitted to the maximum of 
the critical spectrum. The horizontal lines in the lower panel 
exhibit the static susceptibility x'fa ~ * 0) = 1 — / for the 
states with labels n — 2,3, 4, and c, respectively. 



n 3 and then for a frequency interval that is consider- 
ably smaller than the time interval for which Eq. (|2^) de- 
scribes the correlators in Fig. ^. The leading corrections 
in Eq. ( |30| ) are much more important for an adequate 
description of the spectra than for the approximation of 
the correlators. The second correction, Eq. ([nj), is nec- 
essary to describe the plateau for the n = 2 state within 
a 5% error margin. It is also necessary to describe the 
crossover of the spectrum from the plateau towards the 
critical one for n = 3, 4. 

To understand the range of validity of Eq. ( p2[ ) and its 
Fourier transform, one has to compare it with Eq. ( |30| ) 
and its Fourier transform, respectively. This amounts to 
comparing polynomials in ln(f/r) and ln(z/wr). Let us 
restrict ourselves to the dominant terms for the model 
to grasp the essence. Then one can write ((/>(t) — 
f c )/(l - r) = -B\n{t/T)[l - (B 4 /B)\n 3 (t/ T )}. Thus, 
within the error margin S, the leading linear- in-ln(t/r) 
approximation holds for |ln(i/V)| ^ yS{B / B 4 ). For 



the spectrum one gets from Eq. (Al) in leading order 
X »/(l-/ c ) = B(7r/2)[l-4(B 4 /5)ln 3 (l/W)]. Hence, 
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the spectrum is at the plateau within a deviation 5 for 
|ln(l/W)| < $/5(B/ABi). As a result, the range of ap- 
plicability on a logarithmic axis shrinks by a factor 
if one transforms from the time domain to the frequency 
domain. 

Because of Eq. Q, S[ln n (t/r)] {z) = ln"(i/zr) - 
n r yln n ~ 1 (i/zT)+- ■ ■ . In leading approximation for t — * oo 
and z — > 0, one finds x'( w ) oc 1 — (f)(l/u>) whenever </>(£) 
is a polynomial in ln(t/r). This explains the lower panel 
of Fig. |] as a different representation of Figs. ^ and In 
particular, the linear- ln(tj) parts in Fig. |] are of a similar 
size as the linear- ln(t/r) parts in Fig. H. 
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FIG. 5: Correlators <j>(t) for the one-component model de- 
fined in Sec. |iy] for states located on the line ei = —0.0182. 
The states with labels n' = 1-4 have the coordinates 
(vi,va) = (1.1169,2.7141), (1.0669,3.1641), (1.0169,3.6141), 
(0.9669, 4.0641), respectively, and they are marked by tri- 
angles in Fig. hi The state labeled nf — 2 is identical 
with the state discussed in Figs. |l|, ^, ^ with label n — 
2. The states with labels 5, 6, and 7 have the coordi- 
nates (0.9599,4.1271), (0.9569,4.1541), and (0.9549,4.1721), 
respectively. The straight line through the states 1-7 
crosses the liquid-glass-transition curve at the state V^ c = 
(0.95466, 4.17407), where the critical glass-form factor has the 
value f^ c — 0.520. The exponent parameter for the ^-glass- 
transition singularity is A = 0.719 implying a critical expo- 
nent a = 0.318 and a von Schweidler exponent b = 0.608. The 
critical decay law <f>(t) — /' 2 - )c oc t~ a and von Schweidler's law 
4>(t) — f^ c oc — t b are shown by dotted lines labeled a and 
b, respectively; the constants of proportionality are fitted to 
curve 7. The dashed-dotted curve extends the von Schweidler 
expansion for curve 7 to 4>(t) = / (2)c - {t/rf + 1.48(i/f ) 26 . 
The horizontal lines mark the critical glass-form factors f^ c 
and / c , respectively. The dotted straight lines and the dashed 
curves are the leadi ng a symptotic laws, Eq. (p2j), and the lead- 
ing correction, Eq. (130) , respectively. 



Let us consider the states labeled n' = 1-3 and shown 
by triangles in Fig. [l] in order to analyze the implica- 
tions of the correction term in Eq. ( |30| ) proportional to 
£?2- These states are chosen on the line t\ = —0.0182 
and state n' = 2 is identical with state n — 2 consid- 
ered in Fig. H as example for B-i = 0. Figure |^ exhibits 
the correlators together with their approximations. For 



£>2 > 0, the </>(t)-versus-log(t) diagram is concave for all 
times outside the transient, since a parabola with pos- 
itive curvature is added to the leading linear variation 
described by Eq. (^2|). The formula with the leading 
correction describes the complete structural relaxation, 
except for the very last piece for the approach to the 
long-time limit /, as shown by curve n' = 1. This obser- 
vation also holds for cases with £>2 < as is demonstrated 
for state n' — 3. However, for negative B2, the <j>(t)- 
versus-log(t) curve exhibits two inflection points because 
4>(t) crosses the critical form factor f c with negative cur- 
vature. Since the (/>(i)-versus-log(i) curve is convex for 
(j>(t) w f c , it has to have an inflection point for <p(t) < f c 
in order to approach the exponential, i.e. concave, long- 
time asymptote. It has to exhibit an inflection point also 
for </>(t) > f c in order to approach the concave criti- 
cal correlator for short times. The described alternation 
of convex and concave parts is identical to the behavior 
discussed earlier for the MCT correlators for states near 
an A2 bifurcation J22|, [2(|. But, contrary to the char- 
acteristic decay pattern found for the MCT liquid-glass 
transition, curve n' = 3 does not show a two step relax- 
ation scenario, even though there is a huge stretching of 
the dynamics. For the decay from 0.80 to 0.05 a dynam- 
ical window of 5 orders of magnitude is required. Within 
this large window, the correlator follows closely the law 
4>(t) cx ln(i/V cff ). 

The qualitative features described above for state n' = 

3 are more pronounced for state nf = 4, since Bi is de- 
creased to larger negative values. The relaxation curve 

4 has the form expected for states near a liquid-glass 
transition. To corroborate this statement, further states 
5-7 are considered on the line ei = —0.0182 between the 
state 4 and the intersection V~( 2 ) c of this line with the 
liquid-glass transition curve. The transition point V^ c 
is characterized by a critical glass-form factor f^ c > f c . 
The decay of the correlator from the value f^ c to zero 
is the corresponding a-process. Its initial part is de- 
scribed by von Schweidler's power law, as indicated in 
Fig. I for curve n' — 7 by the dotted line. In this case, 
von Schweidler's law accounts for the decay from f^ c 
to about 0.45, i.e. for about 15% of the a-relaxation. 
The analytical description of the a-process can be ex- 
panded by using the extension of von Schweidler's law, 
(t>{t) - f^ c oc ~{t/f) b + B{t/f) 2b |26|, as shown by the 
dashed-dotted line. Asymptotically, the a-process obeys 
the superposition principle: 4>{t) — cj>(t/T a ), where <j> is 
the control-parameter-independent shape function. The 
c6(t)-versus-log(£) curves for the a-process can be super- 
imposed by rescaling the time, i.e. by shifts parallel to 
the log(i)-axis. The reader can check by herself that the 
curves n! = 4-7 have the same shape for <j){t) < f^ c . 
Outside the transient for <p(t) > /'■ 2 - >c , the correlator 
follows the critical decay law for the fold bifurcation 
4>(t) — f^ c oc l/t a , as is also demonstrated for curve 
7. The results for states 4-7 exemplify the well under- 
stood scenario for the evolution of structural relaxation 
near a liquid-glass transition. The formula (p0|) provides 
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an accurate description of 60% of the a-process. 

Comparison of the results for states n' = 1-3 with 
the second-correction formula based on Eq. (^) yields 
the same conclusions as discussed above in connection 
with Fig. H The second-correction formula does not al- 
ter seriously the fit quality for the long-time part of the 
curves n! = 4-6 in Fig. ||. However, for cf>(t) w f^ c , 
the extended formula yields slightly worse results than 
Eq. (p0|). This is so, since for <p(t) > f^ c , the dynamics 
is governed by the A^ singularity V~( 2 ) c whose existence is 
ignored in the expansions near the higher-order singular- 
ity V°. The number f^ c — f c marks the limit where the 
expansion in the small parameter <j>(t) — f c makes sense. 
The opposite conclusion holds for the description of the 
a-process for <j>(t) w f c . Von Schweidler's law results 
from an expansion for states V near V^ 2 ' c in terms of the 
small parameter f@) c — <fi(t). This number becomes too 
large if <fi(t) as f c . It is the dynamics dominated by the 
higher order glass-transition singularity V c that ruins the 
relevance of the expansion resulting in the von Schweidler 
law. The stretching of the a-process connected with the 
transition of V^ c is larger than estimated by von Schwei- 
dler's law, because of the logarithmic-decay effects. 



scribed well by Eq. (|3C|) and it causes a strong a-process 
stretching. The a-peak width at half of the a-peak hight 
is about 2.5 decades. Von Schweidler's asymptotic law 
is irrelevant for the description of the a-peak for states 
n' = 4-6. 



V. RELAXATION FORMULAS FOR STATES 
NEAR HIGHER-ORDER GLASS-TRANSITION 
SINGULARITIES 

In this section, the generalizations of Eqs. (|2^) and 
( |30| ) shall be derived for the asymptotic expansion of the 
solutions dealing with an arbitrary number M of corre- 
lators. In particular, the general formulas for /12, fJ-3, C 
and for the separation parameters e\{V) and (V) shall 
be o btained. The starting formulas are Eqs. (10a) and 
(10t). The subtlety of the problem is the treatment of 



the singular M by M matrix [5, 



qk 



1« c l 
l qk 



A. Asymptotic expansion of the equations of 
motion 




FIG. 6: Susceptibility spectra x"^) f° r the correlators in 
Fig. [| for the labels n' — 1-6. The dashed lines are obtained 
from the leading-correction formula, Eq. (^(i|). The dotted line 
with label a is the critical spectrum proportional to uj a and 
dotted line with label b is the von Schweidler law %' '(lu) oc w~ h 
for state n! = 6. The dotted horizontal line corresponds to 
the spectrum of the leading approximation, Eqs. (E^, p2[ ), 
x"(w) = (1 — / c )v / — 6ei, that is shared by all states n = 1-6. 



Figure ^| exhibits the susceptibility spectra calculated 
from correlators discussed in Fig. |[ The results for the 
states n' = 1, 2, and 3 exhibit the evolution of an a- 
peak if the states cross the line B2 = in the phase dia- 
gram of Fig. |l|. The leading-correction formula, Eq. (30), 
describes this scenario qualitatively. The spectra for 
states n' — 4, 5, and 6 exhibit the superposition princi- 
ple for the a-peak of the susceptibility spectra. The ^3- 
dynamics causes a high-frequency wing of the a-peak fol- 
lowing closely a linear variation with the logarithm of the 
frequency: x"(w) <x — ln(wr). This phenomenon is de- 



The left and right eigenvectors of matrix A^ c for the 
maximum eigenvalue E c — 1 shall be denoted by a* k and 
dfc, k — 1, . . . , M, respectively. According to the Frobe- 
nius theorems [^3|, one can require a% ^ and dfe ^ 0. 
It will be convenient to fix the eigenvectors uniquely by 



the conditions a 
solubility condition of Eq. (10a) reads 



= 1 and a* a 2 



1 

and its general solution can be written as 
4> q {t) = a q $(t) + 4> q (t) . 



1. The 



(35a) 



(35b) 



The splitting of (f) q (t) in two terms is unique if one im- 
poses the condition J2 q a q 4>q{t) — <?(*)■ The part 4> q {t) 
can be expressed by means of the reduced resolvent R q k 



UK 

qk ' 



S[4> q (t)}{z) = R qk J k (z) 



(35c) 



|(l)c 



It is an elementary task to evaluate from matrix A qk 
the vectors a k , a k and the matrix R qk p3|] . 

Equations ( |l0bj ) and ( j35[) suggest an expansion of <j>(t) 
as Eq. (0) and 



Ut) = Gf\t) + Gf{t) + --- , G q n Ht)=<D(\e\ n / 2 ) , 

(36a) 

• , 4 n Ht) = 0(\e\ n / 2 ) . 

(36b) 
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Here, for example, 



Jf(z) = Af\V)+A% 2 a kl a k2 S[G^\mz 



-a 2 Q S[G (1 \t)} 2 (2 



(37a) 



2{A% 2 a kl a k2 S[G^(t)G^{t)]{z) 
-a 2 q S[G^(t)}(z)S[G^(t)](z)} 

A^(V)a k S[G^(t)](z) 

2{A^ k2 a kl S[G^Ht)G k 2 2 \t)]{z) 
-a q S[G^{t)]{z)S[Gf\t)]{z)} 

Afl c k2k3 a kl a k2 a k3 S[G^ 3 (t)](z) 

-a 3 q S[G^(t)] 3 (z) . (37b) 



The justification of the preceding expansions shall be 
given by demonstrating how the equations can be solved 
recursively. 



B. The leading-order contribution 

The leading-order contribution to the solubility condi- 
tion is obtained by substituting Eq. (37a) into Eq. (35a) 



One arrives at: ei(F)+A<S[G« (t)](z) - S[G^ (t)] 2 (z) 
0. Here A = V a* A 



■•<•.' ^ak^k a ki a k 2 1S the expression for the 
exponent parameter 122, \2& and 



e 1 (V) = J2^ q °Hv) 



(38) 



The z = limit leads to ei(V) + (A - 1)/ (1)2 = 0. Com- 
parison with Eq. (£L5ft yields the conclusion that 



M2 = l-E 



* a( 2 ) c 

a q A q kl k 2 a kl 0-k 2 



(39) 



Equation (41) is the factorization theorem of MCT. In 
leading order, 4> q (t) — f q factorizes in two terms. The fac- 
tor h q is time- and control-parameter-independent and it 
characterizes via its g-dependence the specific correlator. 
The other factor is the function G (1) (i) = ~B\n(t/r). 
This factor is shared by all correlators. It describes 
via B and r the control-parameter dependence and via 
ln(£) the complete time dependence. Within the range 
of validity of Eq. (f4l|), the rescaled correlators 9 (1) = 
[4> q (t) — /f] [ hg are the same for all q. Let us emphasize 
that Eq. ( |fl| ) is an exact limit result for the solutions of 
Eqs.© and (§): 



lim[<?!>„ (ir) 



f c q ]/V^W) = -V6/^h q ln(t) . (42) 



The interval of rescaled times t = t/r, where [4> q (t) — 
fq]/ \/—ci(V) becomes close to the r.h.s of Eq. (^), ex- 
pands beyond any bound if V approaches V c arbitrarily 
close. It will be shown below, how the leading corrections 
for 4> q (t) describe violati ons o f the factor izati on theorem 



Substitution of Eq. (37a) into Eq. (35c) yields the 
leading order contribution to 4> q (t), i.e. the function 



G q 2 \t) in Eq. (|36a|). Equation (IT£i) is used to express 



S[G^\t)} 2 (z) in terms of S[G (1) (t)](z) so that 

G q 2 \t)=X q G^ 2 {t)+Yq{V). 

The amplitude X q is independent of e, 



X„ — R 



qk 



A kk x k 2 a ki a k 2 



Yq(V) = 0(e) and reads 



Y q (V) = R qk 4 0) (^)~ £l (^)4 



C. The leading correction 



(43a) 



(43b) 



(43c) 



This parameter has to be zero according to Eq. ( |T8| ) in 
order for V c to be a higher-oder singularity. For fi2 = 0, 
A = 1, and the equation found for G^ (t) is identical with 
Eq. (fl9|). Thus, ei(V) is the first separation parameter 
and Eqs. (|2(]) and ( pT| ) remain valid. 

Introducing the critical amplitude h q by the same for- 
mula as in the theory for the ^-singularity [^2| 

h q = (1 - r q )a q , (40) 

the leading approximation for the correlators is 

Mt)=fq+hqhBHt/T)} ■ (41) 



If one substitutes Eq. ( 43a ) into Eq. (B7a) , one gets an 



Here, B = -Qe^V) / n = C(|e| 1/2 ). Equation @ de- 
scribes the dynamics up to errors of order c. it is the 
generalization of the logarithmic decay law || to arbi- 
trary MCT models. 



expression for J q ' (z) in terms of the known G^\t) and 
the unknown G^ 2 \t). Therefore, the solubility condition, 
Eq. ( |35a| ) , evaluated up to order e 3 / 2 , yields an equation 
for G^j(t), The latter has the form of Eq. (p3|), where 
also the inhomogeneity is given by Eq. (|2^). This holds 
with the formula 

ei{V) = Y.< A q 1 k\ V W + 2e 1 {V)Y j a q a q X q 
q q 

+ 2 E°; [^qtk^YkAV) - a q Y q {V)] (44) 

q 

for the second separation parameter, and the constants 



C = E a ? [a q X q + a q /2] 



(45) 
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M3 = 2C 



E 



k 2 k 3 a kiak 2 a k3 



2A q 2 ktk 2 a kl X k 2 



(46) 



As a result, Eqs. @ for the function G {2) (t) remain 
valid. 

Combining the results for G^{t), G (2 \t) and G q 2 \t) 



with Eq. ( (36a| ), and this with Eq. (||), one obtains the 
main result of this paper. It describes the correlators up 
to errors of order lei 3 / 2 : 



Using Brownian microscopic dynamics, the equations of 
motion ([ftj) and (§ read for tf> q (t), q = 1, 2: 

T q d t <t> q (t) + 4> q (t) + f m q (t - t')d t 4g(t')dt' =(8Qa) 



(50b) 
(50c) 



mi(i) = v l( j)l{t) +v 2 (f>l(t) 
m 2 (t) = v 5 <t>i(i)4>2{t) . 



<f>S) = (f c q + h) + h q [(-B + B,) In(t/r) 
+ {B 2 + K q B 2 )\n 2 {t/T) 



Here 



The three coupling constants v n ^ shall be considered 
as the components of the control-parameter vector V = 

(V!,V 2 ,V 3 ). 

Let us note convenient equations for the discussion of 
the phase diagram [22j, restricting ourselves to v% > 4. 
Equation (^|) for the second form factor implies f 2 = 
+B 3 In (*/ t) + B 4 In (t/ t)\ (.47) _ l]/^/^, and this result can be used to eliminate 

f 2 in the following expressions. Thus, Eq. (j^) for the 
first form factor, /i/(l — /i) = Vxf 2 + v 2 f 2 , is a linear 
equation for (vi,v 2 ) with coefhcients that are nonlinear 
in fi and v 3 . The same statement holds for Eq. (fijl) for 



/, = (!- /„ C )F 9 



(48) 



is a renormalization of t he gl ass-form factor of order e to 
be calculated from Eq. (43c). The critical amplitu de h 
is defined by Eq. (fib] ) . The parameter B\ from Eq. ( 29t ) 
is a renormalization of order e of the prefactor of the 
logarithmic decay law. The three terms proportional to 
B 2l S3, and -B4, respectively, describe leading deviations 
from the logarithmic decay. They are of order e and fol- 
low from Eqs. ( |29c| ) and (29d). The relative size of these 
deviations is the same for all correlators. This means, 
that these terms imply a modification of the factorization 
theorem, (j) q (t) — (/£ + f q ) = h q G(t), in the sense that 
G(t) = -B\n(t/r) in Eq. Q) is to be generalized by 
the bracket on the r.h.s of Eq. (|30"|). It is solely the con- 
tribution proportional to B 2 = —dei/ir 2 = 0(e) which 
describes a violation of the factorization theorem. It en- 
ters with the correction amplitude 



a singularity which is equivalent to /} /(l — fi ) 2 



■2v. 



(2)c z ,(2)c. 



f 2 z>< ~ (I — f 2 2 ^ c )- These equations can 



be used to express i>j 2 ^ c and v^'^ in terms of v. 
/ : . To case the notation, variables x and y shall be 
introduced as 



(2)c 



,(2)c 



and 



,(2)c 



fi 



(2)c 



2/ • 



One gets 



,(2)c _ 



,(2)c 



3- (2 + x)y 



2(1 



2/) 2 2/(2- 
x 2 y{y 2 



xy) ' 
-2y 3 ) 



2(1 -y) 2 (x 2 y 2 -3XJ/ + 2) 



(51a) 

(51b) 
(51c) 



K q = X q /a q 



(49) 



Its q-dependence expresses the fact that the size of the 
leading corrections depends on the chosen correlator. 
Thus, the range of validity of the universal Eq. (^) is 
not universal. The correction amplitude is to be calcu- 
lated from Eq. (43b). The formula for K q is the same as 
discussed in the theory for the A 2 singularity J2fj. 



VI. RESULTS FOR A TWO-COMPONENT 
MODEL 

The simplest example exhibiting a generic swallow-tail 
bifurcation is given by an M = 2 model with the mode- 
coupling functionals T\\V, f k ] = v 1 f 2 +v 2 f 2 2 7 T 2 [V, f k ] = 
vsfif 2 . This model was motivated originally as trunca- 
tion of microscopic equations of motion for a symmetric 
molten salt |2?| . The model shall be used here in order to 
demonstrate implications of our theory which could not 
be demonstrated for the M — 1 model studied in Sec. IV. 



These equations define the surface of bifurcation singu- 
larities of Eq. (0) in the 3-dimensional parameter space. 
The variables x and y with 4 < x and 1/2 ^ y ^ 3/(2+x) 
serve as surface parameters. The exponent parameter 
A = 1 — fi 2 is determined by 



M2 



(3a: 2 + 6x)y 3 - {x 2 + I8x + 8)y 2 + (6x + 18)y - 6 



(2x 2 + 4x)y 3 - I2xy 2 + (2x + 4)y 



(51d) 



The maximum theorem, mentioned above in connection 
with Eq. (0) , has to be used to identify among the points 

(wj 2 ^, v 2 2 ^ c , v^°) those that are glass-transition singular- 
ities. 

Figure exhibits three cuts through the parameter 
space. The cut i>3 = 20 is typical for sufficiently small val- 
ues of x. The cut through the bifurcation surface yields 
a smooth curve of A 2 glass-transition singularities. The 
bifurcation surface for such deals solely with the 

generic scenario for liquid-glass transitions. 

The cut shown for w 3 = 45 is representative for suffi- 
ciently large values of x. In this case, the cubic numera- 
tor polynomial in Eq. (51d) has two zeros yi(x) < y 2 (x) 
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FIG. 7: Phase diagram for the two-component model denned 
in Sec. [Vl]. The full lines are cuts through the bifurcation 
surface for v$ = 20, «3 = and V3 — 45, respectively. The 
value U3 = 24.7 . . . denotes a coordinate of the A4 singular- 
ity indicated by a star. For 113 > 113, there occur A3 glass- 
transition singularities as indicated for v$ = 45 by an open 
circle. The transition lines exhibit a crossing point shown as 
open diamond. The dotted lines that join in the cusp sin- 
gularity marked by a shaded circle complete the bifurcation 
diagram of Eq. ([?]) , but they have no relevance for the discus- 
sion of the MCT solutions (see text). The dashed line denotes 
the V3 = 45 cut through the surface of vanishing first sepa- 
ration parameter t\. The dashed-dotted lines are the cuts 
V3 — 45 through the surfaces of vanishing leading correction 
term B 2 (q) = B 2 + K q B 2 in Eq. @, q = 1, 2. The crosses 
labeled n — 1, 2, . . . and triangles labeled n' = 1, 2. 3 mark 
states whose dynamics is discussed in Figs. |^, and |j, [t(], re- 
spectively. 



above some yo; they can be evaluated elementarily p8[ . 
The transition line consists of several pieces. The first 
one, obtained for 1/2 ^ y > 2/2(2), is shown in Fig. |?j 



as heavy full line. It starts at v^ c — 4, v^''~ = 
and ends at the A3 singularity marked by a circle. The 
second piece describes bifurcations with /12 < for 
2/2(2) > y > yi(x). It connects the mentioned A3 sin- 
gularity with a second A3 singularity of Eq. (^) that 
is shown as a shaded circle. This piece of the line is 
shown dotted. Decreasing y further, one gets a curve with 
fi2 > that joins the second A3 singularity with the point 
v^ c = 0, Ujj = 3/(2 + x). This line exhibits a crossing 
point with the first line piece mentioned above, that is 
shown as diamond. The part between the second A3 sin- 
gularity and the crossing point is shown dotted, and the 
final piece is shown as light full line. The dotted bifur- 
cation lines and the second A3 singularity are excluded 
from the set of glass-transition singularities because of 
the maximum theorem. These items have been added to 
the figure merely in order to allow the reader to recog- 
nize the familiar swallow-tail scenario The crossing 
point organizes three lines of fold singularities. Between 
the A3 singularity and the crossing point, there is a line 
of glass-glass transitions. The continuation of the line to 
the boundary of the admissible parameter range V2 = 



,(2)c 



deals with liquid-glass transitions. Also the third line be- 
tween the crossing point and the parameter-boundary at 
vi = deals with liquid-glass transitions. Both lines are 
characterized by a discontinuous increase of the correla- 
tors' long-time limits from zero to the positive critical 
glass-form factors fq > 0. 

Decreasing x from large values to smaller ones, 
the two cusp values yi(x) and 2/2(2) approach each 
other. The corresponding parameter vectors V c — 
(wj(x), ^2(2), fg(x)) form curves that approach each 
other with decreasing x and join at a certain value x*: 
yi(x*) = 2/2(2*) = y*. The pair (x*,y*) defines the A4 
singularity for the model. The parameters for this sin- 
gularity are obtained, if the derivative of the numerator 
polynomial in Eq. (51d) is zero for ^2 = 0. This leads to 
(x* - 2)0* - 4)(a;*'^30a;* 3 + 136a;* 2 - 168x* + 88) = 0. 
The elementary solution for the zeros of the quartic poly- 
nomial [^8] determines the coordinates of the swallow- 
tail singularity x* = 24.779392 . . . , y* = 0.24266325 .... 
The cut through the transition surface for V3 — x* is 
shown in Fig. |7] as pair of light full lines joining at the 
A4 singularity which is indicated by a star. 
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FIG. 8: Correlators 
fined in Sec. 
on the cut V3 



VI 



i,a (t) for the two-component model de- 
The states labeled n = 1, . . . , 4 are located 
45 through the surface of vanishing dominant 
correction for the first correlator, -62(1) = B2 + K\B 2 = 0. 
The coupling constants are vf — v\ = 2/4™, v% — V2 = 
0.14907/4™ and the states for n — 1, 2, and 3 are shown in 
Fig. ^ by crosses. The full lines are the solution of Eqs. (poja- 
c). The dotted straight lines show the leading approximation, 
Eq. ( ^l[ ) and the dashed ones the leading correction, Eq. (p7[). 
The long horizontal lines show the critical glass-form factors 
/f = 0.312507, /| = 0.92889 and the short horizontal lines 
shown for the states n = 1, 2 denote the renormalized form 
factors fig + fi t 2 according to Eq. (^). Here and in the 
following figures, the model is used with n = Ti = 1. 
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Figure g demonstrates the validity of the factorization 
theorem for states close enough to a cusp singularity V c 
and its violation for states sufficiently away from it. For 
the A3 singularity with v% = 45, the correction ampli- 
tudes calculated from Eq. ([l9|) are quite different for the 
two correlators: K\ = 0.06857, K 2 = -2.049. Therefore, 
the lines for vanishing dominant correction, i.e. the cut of 
the surfaces B 2 (q) = B 2 + K q B 2 = 0, q = 1, 2, with the 
plane ^3 = 45 are quite different as well, as shown by the 
dashed-dotted lines in Fig. 0. The four states discussed 
in Fig. H are chosen on the surface B 2 (l) = 0. Thus, the 
scenario for the evolution of the ln(i/r) law shown for 
4>i (t) is in qualitative agreement with the one discussed in 
Fig. ||. The states with labels n — 3 and 4 are so close to 
the singularity, that the correction term in Eq. (|47j ) pro- 
portional to -82(2) = 0(e) is not important. As a result, 
the rescaled functions (4> q (t) — f q )/h q , q = 1 and 2 agree 
for the states n = 3 and 4, and the same holds for the cor- 
responding approximations. However, for the states with 
labels n — 1 and 2, the negative coefficient B 2 (2) is so 
large that the </>2(i)-versus-log(i) curve does not exhibit 
the straight line obtained for <\>\ (i)-versus-log(t) diagram. 
Rather, the correlator <t> 2 (t) exhibits changes of the cur- 
vature and inflection points as explained above in Fig. |5| 
for state n' = 3. 

Figure || also exemplifies a problem concerning the 
choice of the time scale r. The complete solution of 
Eqs. (H) and (||) is unique up to the choice of a control- 
parameter independent time scale. The nonlinear cou- 
pling of the correlators of different index q requires scale 
universality. However, if a time scale like r is deduced 
from some approximation to the equation of motion, the 
error of the approximation will result in violations of the 
scale universality for the approximate solutions. Con- 
structing the approximate solutions in Fig. || — and also 
in the upper panel of Fig. ^ — the time r was fixed for the 
leading approximation from 4>\ (t) = and for the lead- 
ing correction from </>i(t) = /f + The errors explained 
lead to offsets for the second correlator: <f> 2 (r) 7^ /| and 
02 (t) 7^ f 2 + Hi respectively, for the two approximations 
studied. This explains, e.g., that the dashed line for 4> 2 {t) 
for the state n = 1 does not coincide with the full one. 
One could also choose r differently, e.g. by requesting 
</> 2 (t) = f 2 + h as done in the lower panel of Fig. ^ 

The transition line which is shown in Fig. [?] by the 
light full and almost horizontal curve for the cut V3 = 45 
intersects the line B 2 (l ) = at some glass-transition sin- 
gularity V^ c = (2.94 . . . , 0.130 . . . , 45.0). For states on 
the line B 2 (l) — that are close enough to this singular- 
ity, one gets the standard liquid-glass transition scenario, 
e.g., the evolution of a plateau of the g (f)-versus-log(t) 

diagram at the critical glass-form factor f!p c and an a- 
process for the decay below this plateau. The universal 
bifurcation results for an A4 singularity require that the 
plateau values are below the critical form factors of the 
nearby A3 singularity: f q c < f q . For the example un- 
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Sec. 
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Correlators for the two-component model defined in 



der discussion, one gets f{ 



(2)<: 



0.0747, /, 



(2)c 



0.7027 



The states with labels n = 1, 2, and 3 are lo- 
cated on the line defined by 113 = 45, #2(1) = and have 
coordinates (vi,v 2 ) = (2.7799,0.1183), (2.9254,0.1292), and 
(2.9391, 0.1302), respectively. They are indicated in Fig. frlby 
triangles and approach the liquid-glass transition point V^ 2 ' c 
with coordinates vf )c = 2.941029, w^ 2)c = 0.130326. State 
n' = 1 is identical to state n — 1 discussed in Fig. H. The 
horizontal lines show the critical glass-form factors fg and 
/q 2 ' c , q = 1, 2, for the cusp singularity V c and the fold sin- 
gularity V^' c , respectively. The liquid-glass transition point 
is connected with an exponent parameter A = 0.603, leading 
to the exponents a — 0.363 and b — 0.807. The critical decay 
laws {4> q {t) — fq 2 ^) = h q 2 ^ c (to/i) a are shown as dotted lines 
labeled a. The von Schweidler laws (4>q(t) — fq 2 ^)/^ 2 ^ oc — t b 
with a time scale fitted for curve n' = 3 are shown as dotted 
lines with label b. The straight dotted line in the upper panel 
exhibits the leading asymptotic law, Eq. (|Il|) for </n ft) and 
state n' = 3; the dashed line shows the result of Eq. (flTf). The 
dashed lines in the lower panel exhibit the leading-correction 
formulas, Eq. (^), for 4>2(t) and states n' = 1 and 3, respec- 
tively. 



and /f = 0.3125, /f = 0.9289. The precursor of the 
liquid-glass transition at V^ c explains the stretched tail 
exhibited in Fig. || for the decay of <f>i(t) below 0.1 for 
the state n = 1. 

To corroborate the discussion of the preceding para- 
graph, the correlators with label n = 1 are reproduced 
as curves with label n' = 1 in Fig. @. Two further 
curves with labels n' = 2 and 3 are added. They re- 
fer to states between state 1 and the transition point 
V^ c as noted in Fig. [j] by triangles. The diagrams for 
<j>i(t) for states 2 and 3 exhibit the two-step- relaxation 
scenario characteristic for an A 2 bifurcation. The decay 
for (f> q (t) < /q 2 ^ c demonstrates the superposition prin- 
ciple for the a-process, and its initial part can be de- 
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scribed by von Schweidler's power law. The decay to- 

(2)c 

wards the plateaus fq for t > 1000 follows the crit- 
ical law for the Ai singularity V^ c . The universal 
laws for the dynamics near a fold bifurcation imply that 
the correlators follow the asymptote of the critical law, 
(4> q (t) — fq 2 ' >c )/h q = (t/t )~ a , for short times down to 
about one decade above the end of the transient dynam- 
ics, i.e. until about t — 10. In particular, for small times, 
the correlator for state n' = 2 should approach the one 
for state n' = 3. However, these features are not exhib- 
ited in Fig. [|. Rather, the t~ a law becomes irrelevant 
for the description of the dynamics below times around 
10 3 , where the critical curve crosses the curves de- 
scribing the logarithmic laws for the A3 singularity. As 
a result, there appears a window between the end of the 
transient and the beginning of the description by the A2- 
singularity results where the correlators are described by 
Eq. (|4^). This window deals with an increase in time 
over about two orders of magnitude. In this window, the 
logarithmic decay processes destroy the manifestation of 
the t~ a law. 

The lower panel of Fig. || demonstrates a further im- 
plication of V c dynamics on the precursors of the liquid- 
glass transition dynamics. Even though the time scale 
for the a-process for states n' = 1 or 2 exceeds the one 
for the transient by factors 10 4 and 10 6 , respectively, the 
correlator </>2 (t) does not exhibit the two-step scenario for 
these states. Rather, there is a large time interval where 
the approach towards the plateau follows closely the 

law (<j)2(t) — fi ) l n (^/ r cff)- This is due to cancella- 
tion of two effects: The asymptotes for the V c dynamics 
and for the V^ c dynamics yield a positive curvature, 
while the onset of the a-process causes a negative one. 
The resulting nearly linear-log(t) variation must not be 
mistaken as the true asymptotic logarithmic law given by 
Eq. @. 

The destruction of the critical decay law of the liquid- 
glass-transition dynamics by the presence of a higher- 
order glass-transition singularity nearby alters the fa- 
mili ar p attern of the susceptibility spectra, as shown in 



Fig. 10 The Debye-peak for the transient dynamics deals 
with the spectra for u> > 0.1, as shown by the peak 
around u> « 1 for x'{ This peak is strongly suppressed 
and shifted to higher frequencies for x'ii^)- There is the 
large frequency regime —4 ^ log 10 w ^ —1, where the 
w a -law is irrelevant for the description of the structural 
relaxation spectrum. Rather, the critical relaxation spec- 
trum of the cusp singularity leads to a high spectral en- 
hancement of x'i ( w ) relative to the io a spectrum; it leads 
to a second structural relaxation peak near u> w 0.01 in 
addition to the low-frequency a-peak. It was discussed 
in connection with states n' = 3 and 4 in Figs. || and ||, 
that the winding of the </>(i)-versus-log(i) curve around 
an effective ln(i)-law for states with B2 < is a precur- 
sor phenomenon of a nearby ^-transition singularity. 
Indeed, the spectrum x'ii 1 ^) exhibits the a-peak of the 
mentioned transition with a maximum for logw ~ —3. 
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FIG. 10: Susceptibility spectra for the correlators shown in 
Fig- | 



Thus, because of i?2(2) < 0, the susceptibility for the 
second correlator exhibits two a-peaks, referring to the 
two parts of the liquid-glass-transition lines discussed in 
Fig. |7[ The low frequency a-peaks shift strongly with 
changes of n' = 1, 2, and 3, since the states are shifted to- 
wards the transition singularity V^ c on one of the lines. 
The high frequency a-peak does not change significantly 
since the distance of the states from the other transition 
line is almost unaltered. As explained in connection with 
Fig. ||, the leading correction formula, Eq. (^7|) , describes 
the high-frequency wing of the second a-peak. 



VII. CONCLUSIONS 

Describing the states of a system by a vector V of con- 
trol parameters, the neighborhood of a glass transition 
singularity V c was characterized by a sequence of sepa- 
ration parameters £i(V), e 2(V): ■ ■ ■ ■ These are smooth 
functions of V that vanish at V c , and they are consid- 
ered as small of order e. The glass-transition singularities 
V c are the bifurcation points of the glass-form factors 
fqiV) = (f> q (t — > 00), i.e. of the long-time limits of the 
correlators (j> q (t). These bifurcations are of the cuspoid 
family Ai, I = 2,3... . The ei, . . . , e;_i are the relevant 
small coefficients specifying the polynomial of degree I 
whose largest zero determines f q (V) for small e. The 
major result is the proof that the solution of the MCT 
equations can be asymptotically expanded in polynomi- 
als P{x) of the logarithm of time, x = \n(t/r), for states 
close to an Ai, I ^ 3, and ei(V) < 0. The leading term of 
order \e\ 1 / 2 yields the ln(</r) law, Eq. ([IT]). The prefac- 
tor for this polynomial of degree 1 is given solely by the 
first separation parameter £1(1^). The leading correc- 
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tion adds a polynomial of degree 4. The coefficients are 
of order |e|, and they are determined by e\{V) and e 2 (V), 
Eqs. (|2^) and (E^). The second correction adds a poly- 
nomial of order 7; the coefficients are of order |e| 3 / 2 and 
determined by £i(V), ea(F), £3(V), and so on. Several 
relaxation scenarios have been identified that are utterly 
different from the MCT scenario for the liquid-glass tran- 
sition. The latter is described by an A2 singularity. 

There are distinguished surfaces in parameter space, 
where the prefactor of the x 2 monomial in the polyno- 
mial P{x) vanishes. For this case, the ln(i/r) law domi- 
nates the dynamics for such times where 4> q (t) « f q (V c ). 
This law may describe the complete decay except for the 
transient and for the final exponential approach of 4> q (t) 
towards its long-time limit, Fig. ^. States near the men- 
tioned surface exhibit slight deformations of the straight 
0(t)-versus-log(i) curve. There is a concave behavior on 
one side of the surface and a winding around the straight 
line with alternating convex and concave parts on the 
other side, as shown for the states n' — 1-3 in Fig. ||. 
The corrections to the leading-order asymptotic results 
depend on the correlator under consideration. The sur- 
faces of dominant ln(i/r) behavior are different for dif- 
ferent correlation functions, as explained in connection 
with Fig. 

Every higher-order glass-transition singularity V c is an 
endpoint of a surface of fold-bifurcation points V^ c with 

fg(V c ) = fq < fq )C = f q {V {2)c )- For states sufficiently 
close to V^ c , one finds the standard transition scenario 
with two-step relaxation described by the interplay of el- 
and /3-scaling laws. The correlators for (j) q (t) ss /° are a 
part of the a-process. Therefore, the logarithmic decay 
laws as formulated by Eq. ([l7]) describe the a-relaxation 
master functions. They reduce the range of validity of 
von Schweidler's power law and cause anomalies of the en- 
relaxation shape functions as shown for the states n' = A 
6 in Figs. || and |[ 

Generically, near a higher-order singularity V c , there 
is a further surface of fold bifurcations that crosses the 
transition surface discussed in the preceding paragraph, 
Fig. 0. As a result, there is the scenario for transition sin- 
gularities of type A2, but now with critical form factors 
smaller than the ones at V c : f!p c < f°. Consequently, 
the logarithmic decay laws are a part of the relaxation 
towards the plateau f q 2 ^ c - They reduce the range of 
applicability of the critical decay and introduce a large 
crossover interval for structural relaxation between the 
end of the transient and the beginning of the transition 
dynamics caused by V^ c , as is demonstrated in Figs. ^ 
and [To] for the correlator <fii(t). In particular, there can 
be a crossover from the transient to a simple ln(t/r) law 
followed by a crossover to a von Schweidler power law. 
This scenario was demonstrated by the experiments re- 
ported in Ref. |ll| and by numerical solutions of MCT 
equations for the square- well liquid in Ref. There 
can also be a susceptibility spectrum for structural re- 
laxation consisting of two peaks, as shown for x'ii 1 ^) m 



Fig. 0. 

The asymptotic expansion also describes the critical 
correlator of the higher-order glass-transition singularity 
outside the transient, Fig. |^. These correlators deal with 
the decay towards /° for control parameters at the sin- 
gularity. For states with e\ > 0, <fr q (t) follows the critical 
decay until close to its intersection with the long-time 
asymptote <p q (t — > 00) — f q > f q - Here it crosses over 
rapidly to the glass- form factor f q . Summarizing, the for- 
mulas of this paper provide a qualitative understanding 
of the decay of the correlations provided the state V of 
the system is close to a higher-order glass-transition sin- 
gularity and the correlator <j> q (t) is close to the glass-form 
factor /° at this singularity. 

Let L t denote the length of the log(t) interval where 
an approximation by one of the polynomials P(ln(t/r)) 
describes the solution for the correlator 4> q (t). Let de- 
note the length of the log(w) interval where the Fourier- 
cosine transform of P(ln(f/r)) leads to a description 
of comparable accuracy for the susceptibility spectrum 
X q (uj). It was explained in connection with Fig. |] that 
L u is considerably smaller than L t . This phenomenon for 
glassy dynamics was discussed earlier for the liquid-glass 
transition p(| , but it is more pronounced for the higher- 
order singularities. This feature for stretched relaxation 
is the reason, why it is more difficult to test asymptotic 
MCT formulas with data for spectra than it is with data 
for correlators in the time domain. 

If there is a higher-order glass-transition singularity 
in a disordered system, there is no generic path for the 
evolution of the structural relaxation. Only a parameter 
surface can be generic for the description of the dynamics 
near a cusp singularity. One has to vary two independent 
parameters to identify an A3 singularity, three parame- 
ters to identify an A4 bifurcation, and so on. We hope 
that the demonstration of all basic scenarios for the dy- 
namics near an A3 singularity will be of use to identify 
such singularities in colloids, if there are any. In such 
case, the derived formulas are elementary enough for data 
fitting. Such fitting might lead to a judgment on the rele- 
vance of the subtle implications of mode-coupling theory 
for the discussion of glass-forming systems. 
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APPENDIX A: LAPLACE TRANSFORMS OF 
LOGARITHMS 

The modification of the Laplace transform, introduced 
in Eq. (||), shall be used to map invertibly functions F(t) 
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of time to functions of the complex frequency z. The 
functions are defined for t > and Im z > 0, respec- 
tively. Euler's second integral for the gamma function 
T(y) implies S[t x }{z) = (i/z) x T(l + x) if as > -1. Differ- 
entiating this identity n times for x = 0, n = 0, 1, 2 . . . , 
one arrives at the formula: 



S[ln n (t)] (z) 



E 

k 



r fc in"- fe (*/z) 



(Al) 



Here (™) = n\/[k\{n - k)\] and T k = d k T(x = l)/dx k . 
One gets in particular Tq = 1 and I\ = —7, where 7 is 
Euler's constant. If ij)(y) denotes the digamma function, 
one can write T'(y) — T(y)ijj(y). Iterating this formula, 
one can express Tk in terms of the first (k— 1) derivatives 
of ip(y) f° r V = 1- The latter are given by the tabulated 
values of the zeta function pq| ; for exa mple, T2 — 
r 2 = C(2) = 7r 2 /6. Implications of Eq. (Al) read with 
n ^ 1, ni ^ 1, n2 ^ 1: 

5[ln"(t)] (*)-S[ln(*)r>) = 



12 



n(n — 1) In' 



n-2 



O + E [r fc -rf]in«-*(iXA2) 



S [ln ni+ ">(i)] (z) - 5 [ln ni (t)] (z)S [ln" 2 (t)] (z) 
(7r 2 /6)nm 2 \n ni+n2 - 2 (i/z) 

ni+n 2 



+ E 

fc=3 



E 



ni 



»2 



I v 

x ln ni+n2 - fc (i/-2) • 



(A3) 



These formulas are needed for the evaluation of the func- 
tion f( 2 \z) in Eq. 



Specializing Eq. (A3) to ri\ = n and n 2 = 1 and using 
the definition of the linear operator T from Eq. (|i|) , one 
gets 

T[m"(t)] (z) = (vr 2 /6)[nln"- 1 ^/z) 

71. 

+ ^ (n - k + l)r„, fc ln<"- fe ) (A4) 

fc=2 

where the coefficients are 

(A5) 

Let us construct polynomials p n (x) o f de gree n = 1, 2, . . . 
obeying Eqs. (|26|). Specializing Eq. (A4) to n = 1 shows 
that one can choose pi(x) = x. Assumi ng t he polynomi- 
als for degree I < n to be known, Eq. (A4) provides the 
formula for the degree n 



Pn{x) = X" - y^r n|fc p w+1 -fcW • 



(A6) 



fc=2 



Thus, the sequence of p n (x) can be constructed re- 
cursively in terms of the coefficients T n ^- To derive 
Eqs. (p9[b-d), one needs: 



2.6160a; + x 2 , (A7a) 
-2.1482a; + 3.9239a; 2 + x 3 , (A7b) 
-12.813a; - 4.2964a; 2 + 5.2319a; 3 + a;lA7c) 



P2(x) 
p 4 (x) 
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